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1. Introduction 

Because water has absorption properties for ion- 
izing radiation similar to biological tissue, it has 
been chosen as the standard reference material for 
radiation therapy [1]. Therefore, the absorbed 
dose' water calorimeter was devised [2] and devel- 
oped [3] for investigative studies concerning its use 
as a reliable standard for the calibration of instru- 
ments. One reason for the seemingly great delay in 
initiating such a calorimeter appears to have been 
concerns about the complexities of convection and 
its effect on the accuracy of measuring absorbed 
dose. Although convection at the point of measure- 
ment can be reduced or eliminated by the applica- 
tion of barriers [3] or by operating the calorimeter 
at 4 °C [4], it is of fundamental interest not only to 
better understand the water calorimeter, but also to 
understand convection in a system where there is 
essentially a free flow of water. The few observa- 
tions reported (in an absorbed dose water calorime- 



1 Absorbed dose is the energy per unit mass imparted to matter 
by ionizing radiation. The SI unit of absorbed dose is the gray 
(symbol Gy). 1 Gy= 1 J/fcg. 



ter) are by the authors of reference [A] (who also 
operated their calorimeter at 30 °C) and by Barnett 
[5]. For a water calorimeter operated at ambient 
temperature or higher, it is valuable to know the 
best operational procedures to employ in order to 
reduce undesirable effects as much as possible, 
"With the help of this knowledge and the use of 
convective barriers in the vicinity of the sensing 
thermistor, convective effects can be minimized to 
the point where they are likely to be insignificant. 
A vast amount of theoretical and experimental 
work on convection has been reported in the field 
of fluid mechanics. The most successfully used in- 
strument for measuring small velocities is the hot- 
wire anemometer (1-3 mm long and 5-10 ftm in 
diameter). However, the calibration of this instru- 
ment is difficult and special apparatus is required. 
The lowest velocity at which it has been calibrated 
was about 8 mm/min [6]. A later wort is also de- 
scribed [7], and several difficulties of its use in wa- 
ter are pointed out [8]. Depending on irradiation 
conditions and position of measurement in a bar- 
rier-free water calorimeter, the convective velocity 
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could vary from zero to about 10 ram/min. Even if 
a small velocity could be measured accurately, 
there would still remain the important problem of 
determining the effects of convection on the accu- 
racy of measurement of absorbed dose with a tem- 
perature-calibrated thermistor. The equilibrium 
temperature rise could be significantly disrupted by 
the motion of the water. 

To circumvent the above difficulties, a thermis- 
tor was used both as a temperature-measuring 
device and as a velocity-measuring device. Its ve- 
locity sensitivity is significantly greater than that of 
the hot-wire anemometer. 

In the present investigation, natural convection 
was simulated by forced convection of water flow- 
ing around a thermistor over a wide range of 
known water velocities and thermistor powers. 
The changes in thermistor temperature were mea- 
sured and are presented as the effect of an equiva- 
lent absorbed dose or absorbed dose rate, quantities 
of direct interest to those in the field of radiation 
dosimetry. A conversion factor can be used to 
rescale the results in terms of temperature. 



2, Sources of Convective Effects 

Dissipation of electrical power in the thermistor 
will cause it to rise to an equilibrium temperature 
above that of the surrounding stagnant water, 
which might be either uniform or nonuniform in 
temperature. Irradiation of the water will produce 
temperature gradients which are the driving force 
of convection, and which may set it into motion. 

Measurement errors due to convection as a re- 
sult of irradiation may then arise in two ways: (1) 
if the background temperature of the water is uni- 
form, the motion of the water may disturb the equi- 
librium temperature rise of the thermistor, and (2) 
if the background temperature of the water is 
nonuniform, an additional disturbance may result 
from a temperature difference between the water 
entering and leaving the region of the thermistor. 
The study of effect (1) is the main purpose of this 
investigation. Applications will be described else- 
where [9], 

Nonuniform temperature distribution in the wa- 
ter may have two causes: (a) a significant variation 
of absorbed dose rate in the vicinity of the thermis- 
tor from previous runs, and (b) conduction of heat 
between the calorimeter and surroundings. Cause 
(b) can be reduced by placing an enclosure around 
the calorimeter to circulate the air and regulate it 
within 0.1 °C. Such a simple enclosure was used 



and has been previously recommended [3]. Any re- 
maining temperature gradients, caused by either (a) 
or (b), can be rapidly and effectively erased by sim- 
ply agitating and circulating the water between 
runs. A significant observed drift is a good indica- 
tion that this needs to be done. 



3. Setup 

Figure 1 shows a sketch of the experimental 
setup. A calibrated thermistor (which formed one 
arm of a Wheatstone bridge) was positioned on the 
vertical axis of a two-piece plastic tube, 3 cm inside 
diameter. The tube was supported in a calorimeter 
container C filled with once-distilled water. Plastic 
tubes connected C to an elevated water reservoir R 
and a column of water where micrometer measure- 
ments were made to its surface. 
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Figure 1, Experimental setup for simulating water convection 
to study the change in thermistor equilibrium temperature and 
the resulting effect on a measurement of absorbed dose or ab- 
sorbed dose rate. The drawing is not to scale. The dimensions 
are in centimeters. 

C was surrounded by thick insulation consisting 
of expanded polystyrene (for further details see 
figs. 3 and 5 in [3]). The desired uniform tempera- 
ture in C was achieved with immersion heaters and 
circulating the water with air, which bubbled up- 
ward along the four vertical edges of C. Thermis- 
tor probes in C, R, and the room sensed 
temperatures, which were displayed with a resolu- 
tion of 0.01 °C. 

The large cross-sectional area of R (1,400 cm 2 ), 
helped to maintain a sufficiently constant flow of 
water. The instantaneous initial and final surface 
velocities 5, in C differed by 1 to 5%, as a result of 
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changes in water level heights (about 3 to 13 mm 
inC). 

Measurements were made of the average surface 
velocities S v , rise or fall. This was then used to 
determine the average axial water velocity A, along 
the tube axis. (This is discussed below in section 6, 
S V =£A V ), S v ranged from 1.3 to 10 mm/min. S, (rise) 
was determined from the initial and final microme- 
ter measurements and a measured time interval, 
from the time when the ball valve V, was suddenly 
opened to when it was closed. S r (fall) was deter- 
mined in a similar manner but by operating ball 
valve V 2 ; and, instead of the micrometer measure- 
ments, it was determined from the measured mass 
of water withdrawn and the cross-sectional areas 
of C (927 cm 1 ) and the water column (16 cm 1 ). The 
water lines contained other valves (not shown), 
which were pre-adjnsted to give desired rates of 
water flow when V, or V 2 was opened. 



4. Preparation 

Initially the water was circulated with air and 
the immersion heaters were turned on until the wa- 
ter was heated to room temperature (regulated 
within about ±0.2 °C). This procedure minimized 
subsequent internal temperature drifts. The temper- 
ature of the stagnant water within the tube lagged 
behind that of the surrounding agitated water. To 
speed approach to uniformity of temperature of the 
entire water volume and of the two-piece tube, the 
upper part of the tube was detached (a few cen- 
timeters above the thermistor plane) and placed in 
a corner of C where water flowed vigorously 
within and without that part of the tube. At the 
same time, the motion of the water reduced the 
time to achieve uniformity of temperature in the 
lower part of the tube (still in position). About 30 
minutes later the two parts were rejoined. Several 
minutes later the air supply was turned off and the 
water became stagnant. During the above proce- 
dure, the agitated water in R had been adjusted in 
temperature until it was observed to be the same as 
that in C. 



5. Procedure 

To assure repeatability of water velocities, the 
heights of the water levels in C and R were re- 
stored prior to a run. For example, to simulate up- 
ward convection, V t was opened, which caused 



water levels in C and R to rise and fall, respec- 
tively. After the run, the amount of water that had 
entered C was withdrawn through V 2 and replaced 
in R. To simulate downward convection, V 2 was 
opened, lowering the water level in C. The water 
withdrawn was then poured into R, and Y\ was 
opened until the water levels were restored. 

A procedure was also followed to assure temper- 
ature stability as much as possible. If successive 
runs were made, following the above procedure 
alone, increased temperature instabilities would re- 
sult because of small temperature differences of the 
water flowing into C. This would result in (1) tem- 
perature drifts caused by heat conducted toward or 
away from the tube axis, and (2) a temperature gra- 
dient along the tube axis, which is of more concern. 
Therefore, when the levels were restored after 
each run, the water in C was agitated for several 
minutes. This assured an essentially uniform tem- 
perature distribution along the outside of the tube 
and an eventually more uniform temperature along 
the tube axis. The motion of the water, however, 
did cause a significant disturbance in the thermistor 
equilibrium temperature, but after agitation equi- 
librium was restored within 10 minutes. 

It is, however, still desirable to reduce the rela- 
tively uniform rate of heat conduction to or away 
from the entire tube axis. The direction of the ob- 
served drift (cooling or heating) indicated that R 
needed a small amount of warm or cool water, fol- 
lowed by agitation; or an eyedropper of hot or ice 
water was directly injected into C and agitated. 
The result was a sluggish and delayed reaction, but 
nevertheless this procedure did significantly help 
to provide thermal stability. The remaining small 
drift signals were quickly reduced by making an 
adjustment with a resistance-capacitance circuit 
across the Wheatstone bridge [10]. As with other 
calorimeters, this circuit proved effective in pro- 
ducing an essentially zero initial drift signal. 



6, Axial Water Velocity 

In figure 1, the water surfaces (inside and outside 
the tube) rise or fall at the same rate, S,. Although 
S v will slowly change (because of the changing 
heights with respect to Vi and V^), we consider the 
ideal situation where S„ is constant. 

During the slow laminar flow of a fluid in a tube, 
the velocity at the wall remains zero. As time in- 
creases from f=0, the velocity increases toward 
and along the tube axis, where it is greatest. At 
positions along the axis (sufficiently far from 
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"end effects" of the water surface and the bottom. 
end of the tube), the instantaneous axial velocity A„ 
will increase from &, to 25,. The vertical velocities 
within the tube build up to a parabolic distribution, 
as a function of radius. A detailed theoretical anal- 
ysis is given by Szymanski (i J J. A figure and a de- 
tailed table are presented in reference [11] which 
give the buildup velocity profiles as. a Function of 
time and radius. The figure is also- presented in a 
more recent reference [12], p. 129. Reference [13], 
among many others, gives additional detailed infor- 
mation on fluid Dow. Laminar flow occurs when 
the Reynolds number 



DV 

ifc?= — <2000 

v 



(0 



In equation. (1} U= inside diameter, V= average 
velocity within the tube (V=S t ) r S„=surface ve- 
locity, and v=kinematic viscosity. For the present 
experiment, D=3 em, -3, = 1 cm/min (the maxi- 
mum velocity used), and v—0.59 crnVrrrin (at 
21 °C, the operating temperature of the calorime- 
ter). Substituting these values in equation (1) gives 
Re =s 5 f well within the range for laminar flow. 

In general, the velocity on the axis tor any size 
tube varies approximately linearly with time from 
the initial value S T up to l.5S„. For the 3-cm- 
cbameler tube used in this experiment, the theoreti- 
cal times required For the buildup to reach 50 and 
■95% of the maximum are 0,53 and 2.1 sain, respec- 
tively. Therefore, if at (=0 the surface velocity is 
S r (and remains constant), then at those respective 
times the velocities on the central axis are 1.50iS v 
and 1.95S„. In reality (as mentioned previously) 5, 
is measured, In the data analysis, S r is treated as an 
assumed constant S„„ Therefore, only S„ will be 
mentioned in the remainder of this paper. 

The equilibrium temperature of the thermistor 
was, therefore, disturbed by a steady increase in 
axial water velocity. One phase of the measure- 
ments required measuring the average rate of tem- 
perature change in the thermistor (which can be 
considered as the effect of an equivalent absorbed 
dose rate) during a measured time interval. This 
was related 1o A, during that interval. A t was com- 
puted by use of S T and an empirical fit to the theo- 
retical data mentioned in reference [1 1], integrated 
over the time interval. 



7. The Thermistor 

The thermistor indicated in figure 1 is a commer- 
cially available, electrically insulated thermistor 
(uninsulated thermistors were extremely noisy). 
The nominal size of the bead was 0.25 mm. It had 
two platinum leads (18 ftm in diameter, 1 mm long) 
soldered to paired nickel alloy leads (76 u,m in di- 
ameter). The bead and the JS-pun diameter wires 
had a thin coating of resin, of variable thickness 
<25 ^m. An improvement was mads by applying 
a thin coat of silicone rubber which significantly 
increased the electrical resistance between those 
components and the surrounding water. The result- 
ing outside diameter of the silicone coating around 
the bead was nearly 0,35 mm. The diameter of the 
commercial insulating sleeve around the paired 76- 
iim wires was 0.23 mm. 

The performance of the thermistor was good 
throughout the experiment. The recordings (shown 
in figure 2 and subsequent figures) indicate good 
signal-to-noise ratios. After each day of measure- 
ment, the water in container C was drained and the 
support tube and its thermistor were removed and 
dried in the air-conditioned room, 

P = 19 (Htf 




CDdiTv 



I, Open 

A., = 5.1 jnm/min 



\l, Clased 



Tigpre 2. Typical response -of a che-rmLsror as a result at upward 

cor-'Je-CtLoiL. 

At a calorimeter operating temperature of 21 C C, 
the thermistor resistance was 2.6 kfl. and its mea- 
sured resistance sensitivity was 3.7%/°C. From the 
measurements of its resistance as a function of tem- 
perature and power, che equilibrium temperature 
rise (A7 1 ) of the thermistor above that of the sur- 
rounding stagnant water was: 



AT*- Ul J> mK 



0) 



where P— power in microwatts. 



606 



Volume 93, Number 5, September-October 19SS 

Journal of Research of the National Bureau of Standards 



The value of AT 1 could be significantly different 
for a different bead size and details of insulation. 
Measured results of convective effects are expected 
to vary directly with AT, Therefore, an experi- 
menter should have a similar knowledge of a ther- 
mistor being used, in order to apply the results of 
this investigation. In general, it is good practice for 
experimenters to give an expression similar to 
eq (2) along with the thermistor powers which 
were used. Otherwise, the meaning of some re- 
ported results may be unconvincing to the readers. 

The temperature calibration of the thermistor 
was also determined from the measurements of its 
resistance as a function of its temperature and 
power. This permitted calorimeter responses to be 
converted to temperature changes, described pre- 
viously in detail [3], 



8. Results 

8.1 Upward Convection 

Figure 2 shows a typical general response as a 
result of forced upward convection. In all cases, it 
was possible to preset the power to within an un- 
certainty of 0.1% and to produce an essentially 
zero initial drift. Valve V, was then suddenly 
opened, causing S r =A r =5.l mrn/min. Convective 
cooling of the thermistor was instantaneous, be- 
cause of its increased heat-transfer coefficient ([12], 
p. 409). Its cooling rate depended on A v and power 
P, 10 u.W for the run shown in figure 2. After 2 
minutes, for this case, the thermistor response ap- 
pears to be near equilibrium. At this time, the theo- 
retical velocity buildup reached 94% of the 
maximum [11], giving a value of A v = 1.94St 
= 9.9 mm/min. Then Vi was closed (rapidly caus- 
ing A, to drop to 0), which caused the thermistor to 
rise in temperature to its initial equilibrium value 
represented by the base line LL'. This is an indica- 
tion that the background temperature of the water 
which left the vicinity of the thermistor was equal 
to that which entered that region. When this is not 
the case, the final trace will undershoot or over- 
shoot LL'. Large effects like this were avoided by 
frequent agitation of the water and the preparation 
procedure (described in sections 1 and 4) to assure 
more uniformity in background temperature. Cor- 
rections were made for any of those small remain- 
ing effects caused by nommiformity. 

The scale on the right in figure 2 indicates the 
measurement sensitivity. The indicated absorbed 
dose or temperature change would have caused a 
0. 1-u.V potential change in the output of the 



Wheatstone bridge, which was powered by a 
1.35- volt mercury cell. The noise shown is approxi- 
mately the greatest encountered for all values of 
power. 

The thermistor response during water flow 
showed a wide range of characteristics, which 
varied with power and water velocity. Figure 2 
shows an initial slope S] followed by a steeper 
slope S 2 . Particularly at large signals caused by 
high powers (when the gain of the amplifier was 
reduced, resulting in little or no observable noise 
fluctuations), the transitional point from Si to S 2 in 
many cases appeared to be clear and abrupt; under 
some conditions there was more than one point of 
inflexion. Also, in many cases S,>S> The 
time to reach apparent equilibrium varied signifi- 
cantly. . 

The somewhat complicated response indicated 
by the recordings required some compromises in 
presenting the measured results. The recordings 
show two main characteristics: (1) a rate of change 
of response, and (2) a maximum deflection d aa *, 
represented by the indicated equilibrium level. At 
d msi . the value of A, assigned was 25,, The quotient 
of drnzJl divided by the corresponding time inter- 
val is a measure of the rate of response up to 50% 
of d mta . If d^n is expressed in terms of equivalent 
dose, then that quotient has the dimensions of a 
dose rate. While strictly speaking it is not an "aver- 
age," it is in fact close to a true average, and it will 
be convenient here to refer to it as the average dose 
rate during the time interval rso%- That time inter- 
val appeared to be a reasonable choice, because be- 
yond that point the signal began to decrease 
rapidly. The equivalent "negative" average ab- 
sorbed dose rates D~ are plotted as a function of A* 
during the time interval f i0% , 

A 1 mK rise in temperature is produced by an 
absorbed dose of 4.18 Gy, neglecting a heat defect. 
This conversion factor can be used to rescale the 
plotted results in terms of temperature. 

Figure 3 shows a plot of individual points deter- 
mined during a daily set of measurements, for 
I 3 =25 fiW. The average D~, from to 50% of 
equilibrium, is plotted as a function of the average 
axial water velocity. (The plots in the figures that 
follow do not show the individual points of mea- 
surements, for the sake of clarity.) Figure 4 shows 
plots of the approximate times for the thermistor to 
reach 50% of its temperature drop to equilibrium 
as a function of the equilibrium axial water 
velocity 2S„ 
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Figure 3. Plot of a daily set of measurements. The average 
"negative" absorbed dose rate t)~, from to 50% of equi- 
librium, is platted as a function of the average axial water veloc- 
ity A, during the interval, *»«. 
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Summaries of the measurements (for P = 5 to 100 
p,W) are shown in figures 5, 6, and 1, An expanded 
scale of figure 5 (from A y = to 4 mm/min) is 
shown in figure 6. At low values of P and A r , t>~ is 
very small. Therefore, the curves in those regions 
should serve only as a guide, because some subjec- 
tivity was used in extending the curves to zero A., 
and t>~. Figure 7 shows plots of the equivalent 
"negative" absorbed dose D~ at equilibrium as a 
function of the equilibrium water velocity 2S r 
along the tube axis. 

The smallest values of surface velocity and 
power were 1.3 mm/min and 5 jliW, respectively. 
Six 2-minute runs were made in this condition. 
There appeared to be no detectable response to 
about 1.5 minutes. At the end of the run (A v = 2.5 
mm/min), D~ was perhaps about 0.04 Gy. This 
was barely detectable, being just above the noise 
level. Figures 5 and 6 can be used in practice to 
estimate roughly the convective velocity produced 
by irradiation, by use of the initial and final drift 
rates of a recorded calorimeter run [9]. 
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Figure 5. Summary of the measurements of the average "nega- 
tive" absorbed dose rate t>~ as a function of the average axial 
water velocity A during the time interval tso%. 
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Figure 4. Plots of the time feoss. for Jhe thermistor to reach 50% 
of its drop to equilibrium as a function af the equilibrium water 
velocity IS, along the tube axis. 
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8.2 Downward Convection 
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The investigation of downward convection 
started with initial values of 5 V =3.5 mm/min and 
P=30 u-W. The convective cooling of the thermis- 
tor had the same general shape as that shown for 
upward convection in figure 2. The same was 
noted when P=40 u,W. Most of the time this ap- 
peared to be the case for P = 50 u,W, but occasion- 
ally (after V 2 was opened) there appeared to be a 
short and barely detectable heating response before 
convective cooling began. At P — 60 nW the con- 
vective heating was clearly visible on nearly all oc- 
casions. It always appeared at higher values of 
power, and the response and duration increased 
with power. The power was varied up to 150 fiW. 

Figures 8, 9, and 10 show the development of the 
recordings when the initial velocity was varied, 
P = 80 u.W. Figure 10 may be compared with fig- 
ure 11, where P = \00 uW. 



Figure 6. Expanded scale, fromJ3„=0 to 4 mm/min, of the plot 
shown in figure 5. 
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Figure 7. Summary of measurements of the equivalent "nega- 
tive" absorbed dose D~ at equilibrium as a function of the equi- 
librium water velocity 2S, along the tube axis. 
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Figure 8. Thermistor response to downward convection, 
P = S0 |nW and A„ = l.& mm/min at t =0. 
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Figure 9, Thermistor response to downward convection, 
P = 80 )iW and A r = 3.$ mm/min at t =0. 
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Figure 10. Thermistor response to downward convection, 
P = 80 u,W and A,=7.0 mm/rain at / =0. 



P = 100 ii\N 




Closed 



Figure 11. Thermistor response to downward convection, 
P= 100 jiW and A, =7.0 mm/min at / =0. 



Figures 12 and 13 show, respectively, summary 
measurements of the time interval t\ and the equiv- 
alent absorbed dose D at the peak of response at 
the end of t h (see fig. 8) as a function of power. 
Figures 12 and 13 can be used for estimating the 
average equivalent absorbed dose rates D during ty. 
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Figure 12. Summary of the measurements of the time t[ to the 
peak of the temperature rise as a function of thermistor power. 
The surface velocity S, is in mm/min. 




100 
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Figure 13. Summary of measurements of the equivalent ab- 
sorbed dose D at the peak, at the end of time t u as a function of 
thermistor power. The surface velocity S-, is in mm/min. 
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A reason for the connective heating is postulated 
by use of the illustrations in figure 14. At low ther- 
mistor powers, figure 14(a), the temperature pat- 
tern around the thermistor is essentia]] y symmetric, 
implying no or little convection. Movement of that 
pattern, either up or down, would increase the net 
temperature gradient in the immediate vicinity 

Temperature Patterns 
Around Thermistor At: 



Low 
Power 



High 
Power 



value of power, the results show a wide variation 
in the apparent plume heights. When S v = 1.7 mm/ 
min, the runs from 100 to 150 uW showed that the 
signal was well above the base line IX', when V 2 
was closed. The results of those runs could not be 
shown, because they do not have the required time 
interval ti+t 2 . When 5, =3.5 mm/min, a large 
change in effect was seen. The change is still signif- 
icant when S v = 7.0 mm/min. A still significant but 
smaller change was noted when S„ = 10,0 mm/min. 
It can be seen that the curves in figure 1 5 do not 
converge rapidly. Therefore, no definite conclu- 
sions can be drawn from this information, which is 
presented merely to show this type of behavior. 





(a) 



(b) 



Figure 14. Illustrations of the equilibrium temperature pat- 
terns around a thermistor. 



of the thermistor to cause it to cool. As the power 
is increased, localized convection will eventually 
occur, giving rise to a plume, figure 14(b). Forced 
upward convection would further increase the net 
temperature gradient, since the bottom of the 
plume would be forced nearer to the thermistor, 
and it would be expected that convective cooling 
of the thermistor would be observed. Downward 
convection will cause the plume to move down rel- 
ative to the thermistor. This will initially result in a 
net decrease in temperature gradient, since the bot- 
tom end of the plume will be at a greater distance 
from the thermistor. This will decrease its heat 
transfer coefficient, causing the thermistor at least 
initially to rise higher in temperature, The magni- 
tude of this effect will depend on the axial velocity, 
At a low velocity a longer time will be required for 
the plume to pass through the thermistor, causing a 
relatively long duration of convective heating. 

By the time the value Vj was closed (figs. S to 
1 1), the plume, figure 14(b), had descended down- 
ward, which caused heating of the water over an 
extended distance below the thermistor. This vol- 
ume of water then ascended upward to cause a 
temperature overshoot before temperature equi- 
librium IX' was attained. 

Figure 15 shows plots of the apparent plume 
height as a function of power. This is the product 
of t^+ti and A v during that interval. For a given 



B 
E 

£ 

'5 

X 

0) 

£ 
E 



0) 

1- 

ra 
o. 
o. 

< 




100 
Thermistor Power, yW 



150 



Figure 15. Summary of measurements of the "apparent" plume 
height as a function of thermistor power. The surface velocity 
jS t is in mm/min. 

The above described test for convection around 
a thermistor as a function of power is much more 
sensitive than that previously described [3], where 
it was concluded that no convection of that type 
was observed. 

83.1, Plume Theory The characteristics of a 
plume from a point heat source have been derived 
theoretically [14]. This is also discussed in [15], pp. 
353-354 and [g], pp. 1 10-115. The basic theoretical 
equations apply to the condition where the rate of 
heat transfer from the source is sufficiently high to 
create a laminar convective plume. The present ex- 
periment indicates that this condition is satisfied 
when the power is at least 50 u,W. 
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The theoretical calculations predict the follow- 
ing for a point source in water and P=100 u,W. 
The velocity along the vertical axis is 7.2 mm/min, 
which remains constant with x, the distance above 
the source. The temperature rise at x = 1 mm is 1.5 
mK. Generally, the velocity and temperature rise 
vary as P ln and Px~ l , respectively. The plume ra- 
dius varies as x W2 P~ w . The rate of water flow 
through a horizontal plane at x = \ mm is 
1.1 g/min; the mass flow rate varies as x, indepen- 
dent of P. This surprising result is due to the fact 
that mass flow is proportional to the product of 
velocity and the square of the radius, which is seen 
to be proportional to x [14]. 



9. Conclusion 

The results of this investigation clearly reveal a 
potential source of inaccuracy in measurement of 
absorbed dose using large thermistor powers, in a 
water calorimeter in which uninhibited convection 
is allowed in the vicinity of the thermistor. Figures 
5 and 6 can be used to estimate attenuation of this 
error by inhibiting the convective velocity with 
barriers and reducing the thermistor power, which 
should perhaps be as low as 2 fiW. The output sen- 
sitivity of a Wheatstone bridge varies as the square 
root of the power, so the sensitivity does not de- 
crease as rapidly as the power. Good signal-to- 
noise ratios can be attained at low power levels, 
particularly with the high dose rates available with 
therapy accelerators. To test the initial choice of 
thermistor power and convective barrier geome- 
try, preliminary irradiation measurements should 
be made at several different power levels. No vari- 
ation in measurements should be observed. 



10. Summary 

A calibrated thermistor was immersed in water. 
Electrical power dissipated in the thermistor 
caused it to rise to an equilibrium temperature 
above its initial temperature. The power was 
varied from 5 to 150 u-W, and the equilibrium tem- 
peratures were disturbed by forced downward and 
upward simulated convection with water velocities 
from 1.3 to 17 mm/min. The velocities were re- 
lated to the rate of temperature change, and also to 
the total temperature change from the original 
equilibrium value to a new equilibrium value, The 
rate of temperature change and temperature 
change were expressed as equivalent absorbed dose 



rate and absorbed dose, either positive or "nega- 
tive." 
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